Transmission through a boundary between monolayer and bilayer graphene 
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The electron transmission between monolayer and bilayer graphene is theoretically studied for 
zigzag and armchair boundaries within an effective-mass scheme. Due to the presence of an evanes- 
cent wave in the bilayer graphene, traveling modes are well connected to each other. The trans- 
mission through the boundary is strongly dependent on the incident angle and the dependence is 
opposite between the K and K' points, leading to valley polarization of transmitted wave. 
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I. INTRODUCTION 

Graphene, the latest addition to the family of two- 
dimensional materials, is distinguished by its unusual 
electron dynamics governed by the Dirac equations- 
Wave functions are characterized by spinor whose orien- 
tation is inextricably linked to the direction of the elec- 
tron momentum in a different manner between monolayer 
and bilayer graphenesi 5 -"— Recently monolayer and bi- 
layer graphenes were fabricated using the method of me- 
chanical exfoliation 6 .^ and epitaxially. 9 ^ The purpose of 
this paper is to study the electron transmission through 
boundary between monolayer and bilayer graphenes and 
show that strong valley polarization is induced in the 
transmission probability through the boundary. 

Transport properties in a monolayer graphene arc 
quite intriguing, and the conductivity with/without a 



magnetic held including the Hall effect 



11.12 



quantum 



corrections to the conductivity^ and the dynamical 
transport^ were theoretically investigated prior to ex- 
periments. The magnetotransport was measured in- 
cluding the integer quantum Hall effect, demonstrat- 
ing the validity of the neutrino description of the elec- 
tronic states. 15 ' 16 Bilayer graphene composed of a pair 
of graphene layer a 6 ' 9 ' 17 ' 18 has a zero-gap structure with 
quadratic dispersion different from a linear dispersion in 
a monolayer graphene Z*±2r2k 

In graphenes, states associated with K and K' points 
or valleys are degenerate. A possible lifting of the 
degeneracy has been experimentally observed in high 
magnetic fields^ and there have been various sugges- 
tions on mechanisms leading to valley splitting and/or 
polarization.— ~— A way to detect valley polarization is 
proposed with the use of a superconducting contact^ 

In a graphene sheet with a finite width, localized edge 
states are formed, when the boundary is in a certain spe- 
cific direction! 35 ' 36 Edge states of monolayer graphene 
ribbons have been a subject of extensive theoretical 
studyi^Zr— The electron transport along the boundary 
has been calculated and characterized by odd number of 
channels in each valley4i When the number of occupied 
subbands is odd, a perfectly conducting channel trans- 



mitting through the ribbon is present^ 5 - - — as in metallic 
carbon nanotubesi 5 - ' 60 ' 61 A way to make valley filtering 
has been proposed with the explicit use of the fact that 
only a single right- and left-going wave can carry current 
at each of the K and K' points.— Recently, edge states 
in bilayer graphene were studie d 63 ' 64 and conductance 
through quantum structures consisting of monolayer and 
bilayer graphenes were calculated! 65 ' 66 

In this paper we study boundary conditions between 
monolayer and bilayer graphenes and calculate the trans- 
mission probability as a function of the electron concen- 
tration and the incident angle of injected wave. In Sec. 
II the treatment of electronic states in a kp scheme is 
briefly reviewed and boundary conditions are derived in 
Sec. III. Valley polarization is shown in Sec. IV under the 
condition that the electron density in both monolayer 
and bilayer regions is the same. Numerical results are 
presented in Sec. V and discussion and short summary 
are given in Sec. VI. Analytic results in the vicinity of the 
Dirac point for zigzag and armchair boundaries are dis- 
cussed in Appendix |A1 and IBl respectively, and the num- 
ber of edges states localized at boundaries is discussed in 
Appendix [Cj 



II. MONOLAYER AND BILAYER GRAPHENE 
A. Monolayer graphene 

Figure [T] shows the structure of graphene, two primitive 
translation vectors a and b, and three vectors (I = 
1,2,3) connecting nearest-neighbor atoms. A unit cell 
contains two carbon atoms denoted by A and B. The 
origin of the coordinates is chosen at a B site, i.e., a 
B site is given by Rb = n a & + n^h and an A site is 
Ra = n a SL+ribh-\-T with n a and rib being integers and r = 
T\ = (a + 2b)/3. In the coordinate system (x',y') fixed 
on the graphene, we have a = a(l, 0), b = a(l/2, v3/2), 
and t = a(0, 1/y/S), where a = 0.246 nm is the lattice 
constant. In the following we start with a tight-binding 
model with a nearest-neighbor hopping integral —70- We 
consider the coordinates (x, y) rotated around the origin 



(a) Zigzag-1 (ZZ1 



(b) Zigzag-2 (ZZ2) 




(c) Armchair-1 (AC1 



(d) Armchair-2 (AC2) 



FIG. 1: (color online) Atomic structure near boundaries between monolayer and bilayer graphene. (a) Zigzag boundaries ZZ1 
(rj = -7r/6) and (b) ZZ2 (77 = tt/6). (c) Armchair boundaries AC1 (r? = 0) and (d) AC2 (7? = tt/3). Red (thick) and green 
(thin) lines represent the top layer with a,2 and B2 sites, and bottom layer with Ai and bi sites, respectively. 



by rj such that the y axis is always along the boundary 
of the bilayer graphene. 

In a monolayer graphene, two bands having approxi- 
mately a linear dispersion cross at corner K and K' points 
of the first Brillouin zone. The wave vectors of the K 
and K' points are given by K = (27r/a)(l/3, l/v3) and 
K' = (27r/a)(2/3, 0), respectively. In a tight-binding 
model, the wave function is written as 



V>(r) 



R=Ra 



^A(R)0(r-R)+ <MR)<Kr-R), (1) 



where </>(r) denotes ap z orbital. The amplitude at atomic 
sites R = Ha or Re satisfies 



eV^(R) 
e-0s(R) = -70 



7o ^2 i/jbCR-ti), 



(=1,2,3 

E 

i=l,2,3 



V^(R + Ti), 



(2) 



where the overlap integral has been neglected for simplic- 
ity. 

For states in the vicinity of the Fermi level e = of 
the graphene, the amplitudes are written as 



^(R) = e lK ' R Ff (R) + e 4 "e iK '- R F^'(R), 
V-s(R) = -ue'ie^if (R) +e iK '- R F|F'(R), 



(3) 



in terms of envelope functions Fjf , F%, F^ , and F^ , 3 
where rj is the angle between the x and x' axes as men- 
tioned before and to = e 27 ™/ 3 . The envelope functions are 
assumed to be slowly-varying in the scale of the lattice 
constant. 

For the K point, the envelope functions satisfy the 
Schrodinger equation: 3 



H F(r) = eF(r), 



(4) 



3 



with 



Ho = 7 



F*(r) 



Pf(r) 
Ff(r) 



(5) 
(6) 



where 7 = v / 3a7o/2 is the band parameter, k± — k x ± 
ifcj,, and k = (k x ,k y ) = — zV is a wave vector operator. 
For states with energy e = s^yk with s = ±1 and k = 



fc^ + fc^, the wave function is given by 



F (r) 



sfc_/fc \ jk r 
1 1 ' 



(7) 



apart from a normalization constant. For the K' point 
the Schrodinger equation is obtained by replacing k y with 
— k y and therefore the wave function by replacing k y with 



B. Bilayer graphene 

We consider a bilayer graphene, which is arranged in 
the AB (Bernal) stacking, as shown in Fig. [TJ A bottom 
layer is denoted as 1 and a top layer denoted as 2. The 
unit cell contains two carbon atoms denoted by Ai and 
bi in layer 1, and a 2 and B2 in layer 2. For the inter- 
layer coupling, we include coupling 71 between vertically 
neighboring atoms bi and a2- As a result, the states as- 
sociated with bi and &2 are pushed away from the Fermi 
level, which is the reason that they are denoted by lower- 
case characters. 

Similar equations of motion can be written down for 
amplitudes at atomic sites with the use of nearest- 
neighbor in-plane hopping integral —70 and inter-layer 
hopping integral 71 . In terms of slowly- varying envelope 
functions, the amplitudes are written as 



V^i(R) 

V>a 2 (R) 



iK R t?K 



J F 1 f 1 (R) + e J "e lK '- R Ff 1 '(R), 
^e iK - R F fe f (R) 



iK' R t?K' 



bl 

iK' R t?K' 



Kx (R), 



F& ( R )> 



^ B2 (R) = w-V^V^Fj^R) + r^e* R Ff 2 (R). 

In the vicinity of the K point, for example, the envelope 
functions satisfy the Schrodinger equatiom 3 i 19 i 23 

n F K (r) = sF K (r), 



(9) 



with 



Al bl a2 

/ 7fc_ 
7&+ 71 




B2 
' 


7fc_ 




(10) 



(11) 



We have two conduction bands (s 
bands (s = — 1) 



±^7i + \/|7i 2 + (7fc) 2 



-1) and valence 



0* = 1.2), (12) 



where the lower and upper signs correspond to j = 1 
and 2, respectively. In the energy range —71 < e < +71, 
in particular, we have a traveling mode corresponding to 

£ s ,i(fc) 



ik 



F K (r) = 



\ 



j(k x - iky) 

7i + N 
j{k x - iky) 

7i + |el 
1 



(13) 



apart from a normalization constant. We have also 
evanescent modes decaying or growing exponentially. 
The wave function of the decaying mode in the positive 
x direction, for example, is given by 



G K ( 



7i - |e| 

j(K X - ky 

7i - kl 

1 



-K x x-\-ik y y 



(14) 



with 



N(7i-N) 



+ kl 



(15) 



For the traveling mode, the four-component vector of 
the wave function for k y < is complex conjugate of that 
for k y > 0. For the evanescent mode, however, the abso- 
lute value of the amplitude is quite asymmetric between 
positive and negative k y . This asymmetry is the ori- 
gin of valley polarization of transmitted wave, as will be 
shown below. Further, the bi and a 2 components of the 
evanescent mode diverge at |e| = 71, showing that the 
B2 component vanishes when being properly renormal- 
ized. This is related to the perfect reflection occurring at 
\s\ = 71 for some boundaries as will be shown below. 

In the vicinity of s = 0, i.e., |e| <Si 71, the Hamiltonian 
can be reduced to a (2,2) form with basis set (Ai, Bi) as 



Ho 



.11 
7i 







P 




(16) 



where functions F a 2 and F^i have been eliminated with 



F k [v) « _:U_F# 2 (r), 
7i 

F*(r) a -l-k+F&iT). 
7i 



(17) 



4 



Corresponding energy eigenvalues are 

-v 2 

e s (k) = s—k 2 . 
K ' 7i 



(18) 



This effective Hamiltonian describes the second-order 
process between Ai and B 2 via a2~bi dimers and re- 
produces the low-energy part of the dispersion quite 
well>£ ' 19 ' 20 ' 23 For the evanescent mode given by Eq. (fT4|) 
with Eq. ()15|) . we can neglect |e| in comparison with 71 
in these equations. 

For the K' point, the Hamiltonian is obtained by the 
replacements k y — > —k y . Therefore, the wavefunctions 
are obtained by changing k y into — k y . 



F^(xa,u) for v = K and K' , The same is applicable to 
Eqs. (2jQ) and (HU), giving (x bl , y) = F*(x bl> y) and 
Fg 2 (xB2,y) = for v — K and K' . Because the en- 
velope functions satisfy first-order differential equations 
(fl~0|) . the boundary conditions are fully specified only by 
their amplitudes at the boundary. Therefore, the slight 
deviation of xa and xbi from x = can safely be ne- 
glected and the boundary conditions between envelope 
functions F^ 1 (r), F h " 1 (r) and Fg 2 (r) in bilayer graphene 
and -F4 (r) and Fg (r) in monolayer graphene are written 
as 



F b \(0,y) =F&(0,y), 

F v B2 (o,y) =0, 



(22) 



III. BOUNDARY CONDITION 

Let us consider a boundary between monolayer and 
bilayer graphene as illustrated in Figs. Q] (a)-(d). The 
boundary is straight in the y direction specified by angle 
77. We have zigzag boundaries in both (a) r\ = —ir/6 
(ZZ1) and (b) ry = 7r/6 (ZZ2), and armchair boundaries 
in both (c) 77 = (AC1) and (d) 77 = tt/3 (AC2). For 
these boundaries, the wave functions of both sides can 
be matched only by those in the vicinity of the K and K' 
points, given by Eqs. ((3|) and (jHJ. In more general cases, 
boundary conditions involve evanescent states away from 
the K and K' points, other than those described by Eqs. 
© and ©, and more elaborate treatment is required to 
derive conditions for the envelope functions.— ~— 



A. Zigzag Boundary: ZZ1 

First, we consider zigzag boundary ZZ1 with r\ = 
— 7r/6, as shown in Fig. [T] (a). For A sites on line 
x = xa = —a/(2y/S), we have condition: 

Vm(Ra) = Vai(Ra), =n(a + b) + r 2 , (19) 

where iPai(Ra) is the wave function extrapolated to R4 
from the bilayer region. For bi sites on line x = x b \ = 0, 
we have 



^s(Rbl) = ipbl(R-bl), Rfci = n(& + b), 



(20) 



where ^B(Rfci) is the wave function extrapolated to R^i 
from the monolayer region. Because of the absence of B2 
sites on line x — xb-2 = —a/V3, we have 



V'B2(Rb2)=0, R B2 =n(a + b) 



T3. 



(21) 



The phase of Bloch functions e lK R,A at the K point 
and e zK Ra at the K' point appearing in Eq. (fTS]) given 
by Eq. ((3]) rapidly oscillates as a function of n with period 
of 3 in a different manner. Therefore, the condition (TT9"|) 
is satisfied if and only if the envelope function of each 
valley is the same along line x = xa, i.e., Fai{xa,v) = 



for v = K and K' . 

The boundary conditions do not cause mixing between 
the K and K' points, leading to the absence of inter-valley 
transmission through the boundary. The transmission 
of electron wave through the boundary can explicitly be 
calculated by considering right- and left-going traveling 
modes (J7J) in the monolayer and traveling modes ([13")) 
and an evanescent mode (|14| decaying in the positive x 
direction in the bilayer. Some of the results are presented 
in Sec. V. 

In order to understand how traveling modes of both 
sides are connected with each other, we consider the en- 
ergy region close to the Dirac point |e| -C 71 in the K 
valley. Envelope functions in bilayer graphene are com- 
posed of traveling waves, to be described by F , and an 
evanescent wave G K . The traveling modes in the bilayer 
side are mainly described by two components F^ and 
Fg 2 , and other components are eliminated by using Eq. 
(|17l) . Because the wave vector in the y direction per- 
pendicular to the boundary is conserved, the wave func- 
tions are written as F^(r) = F^ \x)e lkyV , etc. After the 
evanescent mode given by Eq. (| 14[) being eliminated, we 
have following boundary conditions for traveling modes 



Fg 



(0) 



1 



- 



j(k x - in x ) 



7i 



7 



F2M 

Fg(x) 



(23) 

Details are discussed in Appendix A. The boundary con- 
ditions for the K' point are obtained by replacing k y 



with 



-k y . Note that the conditions now include the first 



derivative of the wave functions in the bilayer side be- 
cause they satisfy second-order differential equation (|16l) . 
In the limit \e\/ji 0, they are reduced to 



0. 



T m 



K 



(0), 



(24) 



The amplitude in the bilayer side is asymmetric with re- 
spect to the sign of k y , i.e., the direction of the incident 
wave, and the asymmetry is opposite between the K and 



5 



K' points. This means that for waves incident on the in- 
terface with oblique angle, transmitted waves have valley 
polarization. 

The second condition of Eq. (|24p . together with Eq. 
([7]), shows that the reflection coefficient becomes tkk ~ 
— 1 and the transmission probability vanishes when an 
electron wave is incident from the monolayer side. On 
the other hand, the first equation of Eq. (|24|) shows that 
the amplitude of the wave transmitted into the bilayer 
side is appreciable unless Fj£ (0) = 0. These somewhat 
contradictory conclusions arise from the fact the trans- 
mission probability is multiplied by the velocity which 
is proportional to k in the bilayer side and therefore is 
much smaller than in the monolayer side. Some examples 
of the wave functions will be shown in Fig. |6] 



B. Zigzag Boundary: ZZ2 

For the zigzag boundary ZZ2 (j] = 7r/6) illustrated in 
Fig. HI (b), boundary conditions become 



V>a 2 (Ra 2 ) -0, 



R,4 = nb + f 2 , 

Rbi = nb, 
R a2 = nb - el, 



(25) 



giving conditions for the envelope functions 

F^(0,y) =FX(0,y), 
F&(0,y) =FH,(p,y), 



(26) 



In the vicinity of the Dirac point |e| <C 71, boundary 
conditions for traveling modes become 



*f (0) 



(l- 







S"f(k x +ik y ) "f(k x ~ik y ) 



7i 7i 
In the limit \e\/-fi — > 0, they are reduced to 




wave functions, we have boundary conditions 



Ff (0) « ( 



1 



(28) 



f#(0) 



0. 



Essential features of the boundary conditions are the 
same as in the case of ZZ1. This fact will be demon- 
strated by approximate but analytical results in Sec. IV 
and by numerical results in Sec. V. 



C. Armchair Boundary 

Next, we consider armchair boundary AC1 (77 = 0) 
shown in Fig. [1] (c). By a proper extrapolation of the 



V>a(Ra) 
Vu(Rai) 
^b(Rs) 
V'b(Rw) 

^B2(RS2) 
V , a2(Ra2) 



= i> A \ (Ra), 

= ^l(R B ), 
= "0&l(Rbl), 

-0, 
= 0, 



R A = n(a+2b) + f 2 , 
Kai = n(a + 2b)+n, 
R B = n(a + 2b) + b, 
R w = n(a+ 2b), 
R B2 = n(a + 2b)-r 3 , 
Ra2 = n(a+2b) + b, 



(29) 



where and Rb are on line x — x 2 = —a/2, Rai 
and Rfci are on x = x\ = 0, and R_b 2 and R q2 are on 
x = x 2 - Because K • (a + 2b) = K' • (a + 2b) = (mod 
27r), the Bloch functions remain constant on lines x = ir 2 
and x = x\. Thus, we have from the first and second 
conditions of Eq. (|2"9")l 



Ff(r) 



F£(r)+F£(r) - * A , 
wF^(r)+F^(r) =«Fj(r) 



Ff (r)| 
+-Ff (r)l = 



(30) 



a;— a;i ' 



respectively. Note that the slight deviation of x± and a; 2 
from x — can safely be neglected from the same ar- 
gument for zigzag boundary. Because envelope functions 
are slowly varying in the scale of a lattice constant, both 
conditions are satisfied, if and only if they are the same 
within each valley. Exactly the same argument is appli- 
cable to the third and fourth conditions. For the fifth 



and sixth conditions of 
and e *(K'-K).R a2 = ^ 



we use e 



;(k'-k)-Rb 2 



Then, the boundary conditions 



for the envelope functions are summarized as 



Fi 2 (0,y)+F^(0,y) 



F a K 2 (0,y) 



F v B {^y), 
0, 

0. 



(31) 



Armchair boundary AC2 of n = 7r/3 is illustrated in 
Fig- [2(d). In a similar manner, the boundary conditions 
are obtained as 



Fl^y) =F%(0,y), 



"F«(0,y) 



F 



F&(o,y) 

K' 1 



a2 (0,2/) 
(0,2/) 



F K ' 



F^(o,y), 
= 0, 
= 0. 



(32) 



These conditions are converted into those of AC1 (|3Tj) 
by changing the relative phases of the envelope functions 
for the K and K' points. Therefore, there is no differ- 
ence between transmission probabilities, etc. of AC1 and 
AC2 within the present k-p scheme, although actual wave 
functions ^(R), etc. may be different. 

Inter-valley mixing occurs at the armchair boundary 
in contrast to the zigzag boundaries. Effective boundary 
conditions in the vicinity of the Dirac point, |e| <C 71, 
can be derived in a manner similar to those for the zigzag 
boundaries and the results are presented in Appendix [B] 
There, we show that the conditions are essentially similar 
except for the presence of small inter-valley mixing. In 
fact, in the limit of k — > 0, an injected wave is perfectly 
reflected within each valley, i.e., tkk — —1 and tk'k = 
for wave incident in the K valley and the transmission 
increases with energy as for the zigzag boundaries. 



K 


K' 






ky > ky < 


ky ky 


< 


Amplitude 


ZZl 1 (1) 


1 (1) 





Ai, a 2 


ZZ2 1 





1 


B 2 


AC1/AC2 










TABLE I: The number of e = edge states present in the 
bilayer graphene localized at the boundaries of monolayer and 
bilayer graphenes. The carbon sites where the wave function 
has nonzero amplitude are shown in the rightmost column. 
The number in parenthesis denotes that of perfectly reflecting 
states present at |e| =71. 



D. Edge States and Perfectly Reflecting States 

As in monolayer and bilayer graphenes, there exists an 
edge state at e = with amplitude only in the bilayer re- 
gion localized at the boundary for zigzag boundaries (ZZl 
and ZZ2) and no edge state for armchair boundaries as 
is shown in Table HI The details on the derivation are 
discussed in Appendix [C] These edge states do not play 
important roles in the transmission through the bound- 
aries because the transmission is possible only away from 
e = 0. In Appendix O further, we show that at |e| = 71 
we have perfect reflection in the region k y < at the 
K point and k v > at the K' point only for boundary 
ZZl. These states are also included in TableQ] This spe- 
cial feature of ZZl clearly appears in numerical results 
presented in Sec. V. 



IV. VALLEY POLARIZATION 

We consider electron transmission between a mono- 
layer and bilayer graphene with same electron concentra- 
tion. This is realized when the electron density is changed 
by a gate voltage. In the presence of electric field due to 
gate, the symmetry between the top and bottom layers 
of a bilayer graphene is broken and a small band gap 
can openi 21 i 71 i 72 ' 75 This small gap will be completely ne- 
glected in the following, because we are interested in the 
essential feature of the transmission property. Besides, 
the Fermi level always lies away from the gap and the 
asymmetry can be controlled by the field due to an extra 
gate. The electron density higher than 7/2/71 > V2 can 
experimentally be achieved by various methods i 73 ' 74 

Electron wave with wave vector k and positive group 
velocity in the x direction at Fermi energy sp is injected 
from the K valley in the monolayer side at the Fermi 
level. For incident angle 9 (—tt/2 < 9 < n/2), we have 



k x = skcos9 and k y = sksin9 with k — |k| for the 
incident wave. The wave is reflected in the direction it— 9. 

When jk/ji < V2, only a single conduction band is 
occupied in the bilayer. In this case, the wave trans- 
mitted into the bilayer has the same wave vector k, i.e., 
there is no refraction. When 7/0/71 > V2, two bands 
are occupied by electrons in the bilayer, giving rise to 




' k-H f -Hsh^ 



(b) 

ykly^>2 m 



n v S Fermi Circle ^ ^' 

Monolayer Bilayer 

FIG. 2: Some examples of the alignment of energy bands 
and corresponding Fermi circles in monolayer and bilayer 
graphene under the condition of equal electron concentration. 

two Fermi circles. In this case, the number of transmit- 
ted waves changes from two to one with the increase of 
9 and the total reflection occurs for sufficiently large 9. 
This is illustrated in Fig. [2J 

In the energy region close to the Dirac point 7/0/71 — > 
0, a simple expression can be obtained for the amplitude 
/ of the transmitted wave for incident wave given by 
Eq. ([7]). The details on the derivation are discussed in 
Appendix A. The result is 



2s cos# 



(sin 61- sy/l + sin 2 9) 2 



(33) 



Because the velocity is j/h in the monolayer and 
2j 2 k/jih in the bilayer, the transmission probability is 
proportional to fc|/| 2 - Therefore, it vanishes for k — in 
agreement with tkk = — 1 as discussed in the previous 
section and increases in proportion to k. Further, it takes 
a maximum at 9 = s9q, with 



1 

7! 



0.196tt. 



(34) 



For the K' point the amplitude is obtained by replacing 
9 with —9. The valley polarization 62 of the transmitted 
wave becomes 



P 



T K -T k > 2 sin 9 y/l + sin 2 9 



(35) 



T K + T K , l + 2sin^6> 

where Tk and Tk' are transmission probability into K 
and K' valley, respectively. The valley polarization in- 
creases with incident angle 9, up to P — ±2\/2~/3 « 
±0.94 at 9 = ±7r/2. 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



Angle (units of it/2) Angle (units of ji/2) Angle (units of nl2) 



FIG. 3: (color online) Calculated transmission probabilities of the K valley as a function of incident angle 6 for several charge 
densities specified by k in the monolayer graphene. (a) A zigzag boundary (ZZ1) with 77 = — 7r/6. (b) A zigzag boundary (ZZ2) 
with 77 = +7r/6. (c) An armchair boundary (AC1 and AC2) with 77 = and +7r/3. The vertical dot-dot-dashed line shows the 
maximum-angle in the limit k — > 0. 



For a ZZ2 boundary, from the first equation of (|28|) . 
the amplitude is calculated as 

2s cos 9 



f 



We have 



_ e -2i8 _|_ j e -i0( sm _ sVl + sin 2 



/ = V2f e 



i(9+a) 



-s tan 



+ sm 



(36) 



(37) 



cos 9 

where fo is denned in Eq. (|33[) for ZZ1 boundary. There- 
fore, a maximum transmission also occurs at 9 = sdo 
for the K point and — s9o for the K' point. For arm- 
chair boundaries, the analytic expression of the ampli- 
tude is presented in Appendix [Bj It shows that the inter- 
valley transmission probability between K and K' is 1/5 
of the intra-valley transmission for perpendicularly inci- 
dent wave (9 = 0) near the Dirac point and that max- 
imum transmission occurs at 9 w s x 0.1797T for the K 
point and 9 ~ — s x 0.179-/T for the K' point. 

The valley polarization completely disappears when 
two traveling waves are involved in the transmission in 
the bilayer, i.e., for small incident angles in the case 
7/0/71 > v2' In this case the wave functions for — k y are 
simply obtained by taking complex conjugate of those for 
k y in both monolayer and bilayer graphenes and there- 
fore the reflection and transmission coefficients for —9 are 
related to those for +9 through complex conjugate. Con- 
sequently, the transmission and reflection probabilities 
become symmetric about 9 = 0, as will be demonstrated 
in the next section. Asymmetry reappears at large inci- 
dent angle for which transmission into a single traveling 
wave is allowed. 



V. NUMERICAL RESULTS 

Figure [3] shows some examples of calculated transmis- 
sion probability as a function of incident angle for (a) 
zigzag boundary ZZ1 with r\ = —tt/6, (b) zigzag ZZ2 
with 77 = +7r/6, and (c) armchair (AC1 and AC2) with 
77 = and +7r/3. The electron density is specified by 
k corresponding to the Fermi energy in the monolayer 
and the results in the low-density regime 7A/71 < V2 
are shown. The transmission probability varies strongly 
as a function of the incident angle and its maximum ap- 
pears at an angle deviating from the vertical direction. 
This asymmetry is opposite between the K and K' points, 
showing that strong valley polarization can be induced 
across the interface of monolayer and bilayer graphenes. 
Except in the high-concentration region 7/0/71 ~ y2, the 
valley polarization is similar for different boundaries. 

Figure [4] shows the total transmission probability in 
the high-density region 7/0/71 > y/2. The transmission 
probability depends strongly on boundaries. In fact, at 
the bottom of the first excited conduction band, i.e., 
^7/71 = v2, it completely vanishes in the region 9 < 
for ZZ1, but not for ZZ2 and armchair boundaries. This 
vanishing transmission at s = 71 for boundary ZZ1 is 
closely related to the presence of a perfectly reflecting 
state in the region k y < and k y > for the K and 
K' point, respectively, as discussed in Sec. III.D and Ap- 
pendix [cj 

This can also be understood directly from the behavior 
of the evanescent mode given by (1 14|) and the boundary 
condition. At |e| ~ 71, the amplitude of the evanescent 
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FIG. 4: (color online) Calculated transmission probabilities for yk/yi > v2, corresponding to Fig. [3] Fermi lines in monolayer 
and bilayer graphene are depicted in the inset. Two transmitted waves are present in the bilayer in the region between two 
upward arrows and no transmission is allowed outside the region denoted by downward arrows. 




0.0 0.5 1.0 1.5 

Fermi Wave Vector: k (units of y-|/y) 



FIG. 5: (color online) Calculated transmission and reflection 
probabilities through armchair boundary AC1 as a function 
of the Fermi wave length in the monolayer graphene. Thin 
lines show analytic results valid for 7/0/71 <C 1. 



mode is nonzero at 02 sites and vanishes at B2 sites, 
and thus it cannot contribute to boundary condition of 
ZZ1, Fb2(0,u) — 0. As a result, the condition should 
be satisfied by the traveling mode alone, leading to the 



vanishing amplitude of the transmitted wave. For ZZ2, 
on the other hand, the condition F a 2(0,y) = is easily 
satisfied even for nonzero amplitude of the traveling mode 
because of the evanescent mode, leading to appreciable 
transmission. 

When 7/0/71 > V2, the first excited conduction band 
crosses the Fermi energy and thus the second traveling 
mode opens for small incident angles between upward 
arrows. In this case, the transmission probability is sym- 
metric about = 0, causing no valley polarization, as 
has been discussed in the previous section. For large in- 
cident angles (outside the downward arrows), there is no 
traveling mode in the bilayer graphene and therefore the 
transmission probability vanishes. 

Figure [5] shows transmission and reflection probabili- 
ties through the armchair boundary. Inter-valley trans- 
mission and reflection probabilities are much smaller than 
the intra-valley probabilities when the electron density 
is sufficiently small, but slowly increase with energy and 
become comparable to intra-valley probabilities when the 
Fermi level reaches the bottom of the first excited con- 
duction band. 

Figure [6] shows some examples of the wave function as 
a function of position for a zigzag boundary ZZ1 with 
7] = — 7r/6. The energy is chosen to be sufficiently small, 
the incident angle = 0. We note that Fb(0) in the 
monolayer graphene becomes vanishingly small and con- 
sequently Fa (0) ~ 2 in agreement with the discussion in 
Sec. IV. Further, the boundary conditions (f22j) are sat- 
isfied by the presence of considerable amplitude of the 
evanescent mode. In fact, the spatially-varying ampli- 
tude in the region x > mostly consists of the evanescent 



0.0 




Bilayer 



Zigzag-1 , r\=-n/6 
6=0 

k=0.01yi/y 



2000 
Position (units of a) 



4000 



FIG. 6: (color online) Calculated envelope function around a 
zigzag boundary ZZ1 (77 = — 7r/6) for wave incident normal 
to the boundary. The amplitude varies as a function of x 
in the monolayer (a; < 0) due to interference of incident and 
reflected waves. 



mode. 



VI. DISCUSSION AND CONCLUSION 

Explicit numerical calculations have been performed 
within the model of uniform charge density on both 
monolayer and bilayer regions. In this model, the energy 
measured from the Dirac point can be slightly different 
between the layers when the electron density becomes 
nonzero (see Fig. [5]). In actual systems, this may be re- 
alized by the presence of small potential variation in the 
vicinity of the boundary, which should be determined in 
a self-consistent manner. The essential features of the 
results that envelope functions are well connected at the 
boundary and that strong valley polarization occurs due 
to the boundary transmission are expected to be inde- 
pendent of the presence of such small perturbations. 

We can also consider the case that the kinetic energy 
of the incident and transmitted waves is the same be- 
tween two regions. This is realized, for example, when a 
hot electron above the Fermi sea is injected. The trans- 
mission is understood in the same manner, but there ap- 
pears some significant difference because of the difference 
in the wave vector of the monolayer and bilayer, in par- 
ticular, when the Fermi level lies in the vicinity of the 
Dirac point. For a given wave vector k in the monolayer, 
for example, the wave vector becomes k' — y/jik/'y 3> k 
in the bilayer. Because the wave-vector component k y 
parallel to the boundary is conserved, this leads to the 
focusing of the transmitted wave into the vertical direc- 



FIG. 7: Schematic illustration of ribbon and prism shaped 
region of bilayer graphene with zigzag boundary. In a ribbon- 
shaped bilayer with zigzag structure, boundaries always con- 
sist of a pair of ZZ1 and ZZ2, but in the case of equilateral 
triangles, all boundaries consist of either ZZ1 or ZZ2. 



tion, i.e., \9'\ < arcsin le/711 1 / 4 , where 6' is the angle of 
the transmitted wave. Further, we have k x w k' 3> k y , 
showing that k y can be neglected in Eqs. (j24j) and (|28|) . 
Then, the transmission is nearly independent of the in- 
cident angle and the valley polarization is considerably 
reduced. 

The valley polarization of waves transmitted through 
a single boundary is reduced when waves go through 
a ribbon-shaped narrow bilayer region sandwiched by 
monolayer graphenes, as shown in Fig. [T^a). The rea- 
son lies in the cancellation at two parallel boundaries. 
The time-reversal symmetry gives the relation that the 
transmission probability incident from the monolayer at 
the K point with angle 9 is the same as that incident from 
the bilayer at the K' point in the reverse direction, i.e., 
|ig3£(0)| 2 = \t%? K ,(9)\ 2 , where 'BM' and 'MB' stand for 
waves transmitted from monolayer to bilayer and from 
bilayer to monolayer, respectively. Let us consider a hy- 
pothetical ribbon consisting only of ZZ1 boundary. With 
the use of the symmetry |i^|(6»)| 2 = \t^\, (-6»)| 2 , the 
total transmission probability through the bilayer ribbon 
is proportional to |i^|:(0)| 2 x \t™{8)\ 2 = \tff,(-9)\ 2 x 
|i^-(0)| 2 , when interference effects are neglected. The 
result is independent of K and K' points. 

Actually, zigzag bilayer ribbons always consist of a pair 
of ZZ1 and ZZ2 boundaries as shown in Fig. [TJa), giv- 
ing different amount of valley polarization. Therefore, 
the cancellation is not complete and certain amount of 
valley polarization remains after transmission through a 
ribbon except in the vicinity of the Dirac point |e| <C 71, 
where the transmission probabilities across ZZ1 and ZZ2 
are different only by factor two, leading to the com- 
plete cancellation. This cancellation is reduced for two 
boundaries not parallel to each other and the polarization 
can be enhanced, for example, when waves go through 
a triangular-shape bilayer island formed in a monolayer 
graphene as shown in Fig.[7][b). 

Boundary conditions for edges of monolayer graphene 
with more general forms were discussed previously and 
boundaries were shown to be classified into either arm- 
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chair or zigzag types^ Similar considerations are likely 
to be possible in the present system. For interfaces 
other than zigzag and armchair, however, the full bound- 
ary conditions require the presence of evanescent modes 
which are not described by states in the vicinity of the 
K and K' points given by Eqs. © and dHJ). 67 " 70 This 
problem is left for a future study. 

In conclusion, boundary conditions between mono- 
layer and bilayer graphene have been obtained within an 
effective-mass scheme based on a tight-binding model. 
Evanescent mode decaying exponentially away from the 
boundary plays an important role and as a result the 
traveling modes are strongly connected to each other be- 
tween the monolayer and bilayer graphenes. The trans- 
mission probability can be quite different between K and 
K' states for waves incident in oblique directions, result- 
ing in significant valley polarization of waves transmitted 
through the boundary. 
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The envelope function in the monolayer side consists 
of incident wave in the direction 9 and reflected wave in 
the direction tt — 9, i.e., 



F§{x) 



(A3) 

with reflection coefficient tkk, where we use kx — 
sk cos 9 and k y = sk sin 9 in Eq. ([7]) . Under the condition 
of equal electron density in the monolayer and bilayer re- 
gions, the transmitted wave is written as 



with amplitude /. 

Upon substitution of Eqs. ([Ml and ([ATI) in Eq. 
and 



(A4) 



K x ~ \J\e\lill 2 + « k\/l + 



sin 



we have 



/ = 2cos6>( 



-2i8 



(A5) 



(A6) 



tkk 



1 - — (s cos 9 - i\fl + sin 2 o) f. (A7) 



7i 



The transmission probability is given by \txK\ 2 = 
2(7/2/71) |/| 2 . This satisfy the unitarity condition \tKn\ 2 
+ Irffif | 2 = 1 up to the lowest order in 7/2/71. 



Appendix A: Low Energy Approximation 

In order to understand boundary properties, the 
boundary condition (|22l) is examined in the low energy 
approximation |e| -C 71 ■ The envelope function in bilayer 
graphene consists of traveling wave F^ and evanescent 
wave G K . Then, Eq. ((55J) becomes 



b A 
F K 

r B 



r bl 



=F§ 2 



- aGf 2 , 



(Al) 



with coefficient a. In the low-energy regime, F^f can 
be replaced by F§ 2 with the use of Eq. (fTT|) and the 
evanescent wave (|14|) is approximated by 



kii 




\ 



kn 



-is 



7i 

7( K a; ky) 



7i 



V 



/ 



Eliminating a in Eq. (IA1[) . we immediately have Eq. 



Appendix B: Armchair Boundary 

An armchair boundary AC1, for example, shall be dis- 
cussed in the vicinity of the Dirac point |e| <C 71. After 
elimination of evanescent modes of the K and K' points, 
boundary condition for traveling modes becomes 



Fa(0) 
F b (0) 



f Faa Tab 

IfflA TBS 



Fai(x) 

F B 2(x) 



(Bl) 



x=0 



with 



^AA — 



^AB = — - 



f^X ~\~ ky 


2k x 


K x ~\~ky i 


k x ky 


K X ky 




^X ky 


K X ky 


Kx 


&X 


K X ~\- ky 


K x H~ ky 



K X ky lk X ~\- ky \ 

Kx~\~ky C 2'Kj X 
^ Kx~\~ky lk X ~\~ky 
K j x ^k^y 2*K> X 



(A2) T BA = g 



r\>x Kx 
K*x I 'ky / C t \ 



' {kx ~\-iky) 



{k x i^y) 
~ {kx iky ) 



-iky 



2. ( k * 

7i V kx + iky 



1 7 KZ — k 



2 71 K 



l-K (1 1 



(B2) 
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and 



etc. 



Fa(x) 

In the limit |e|/7i — > 0, they are reduced to 



Ff(x) 



Fa(0) ^t aa f a1 (x)\ x=0 - 

F B (0) «0. 



T A bF B2 (x)\ x=0 , 



(B4) 



This shows that Tkk = — 1 an d r K'K = for electron 
wave incident from the K valley at k — 0, the same as for 
zigzag boundaries. With the increase of k, the transmis- 
sion increases in proportion to k and its amplitude can 
be estimated using the first equation of (|B4|) . Because of 
the presence of off-diagonal elements in V aa and Tab, 
inter-valley mixing occurs at the armchair boundary in 
proportion to k. After some manipulations, the ampli- 
tude / transmitted into K valley and /' into the K' valley 
become 



/ = "7« 



2iB 



s0 (S - e 2i0 - 2zse i Vl + sin 2 6>) 



(B5) 



a 



s tan 



sin 



sin 



For the K point, the transmission probability is propor- 
tional to k\f\ 2 which takes maximum at 9 « 0.1747T, and 
for the K' point k\f'\ 2 which takes maximum at 9q. Anal- 
ysis of the above equations reveals that inter-valley mix- 
ing is 1/5 of the transmission probability within valley 
for perpendicularly incident wave (9 — 0). The total 
probability is given by the sum of them and maximum 
transmission occurs at 9 sa 0.1797T. 



Appendix C: Edge States 

As in edges of monolayer graphene ) 35 i 36 there can 
be edge states localized along a boundary between the 
monolayer and bilayer graphene. An edge state consists 
of evanescent modes exponentially decaying in the nega- 
tive x direction in the monolayer and those decaying in 
the positive x direction in the bilayer. In the following, 
we shall confine ourselves to the case of vanishing electron 
density in both monolayer and bilayer regions. 

In monolayer graphene occupying half space x < 0, a 
relevant evanescent mode with energy e and k y in the 
range \yky\ > |e| has imaginary wave vector in, with 



1 K = -^( 1 k y y-e 2 , (CI) 
and the wave function Q K e - KX + lk vV f or the K point, with 



(C2) 



where s £ and s y denote the sign of e and k y , respectively. 
The wave function for the K' point is obtained by replac- 



(B3) ing k y with — ky 



In bilayer graphene lying in the region x > 0, we can 
have two evanescent modes with wave vector 



7«j = +\/(l k v) 2 ~£ 2 + Sj|e|7i, 

J \ +1 U = 2), 
and wave function Q, K . e - K i x + lk yV ; with 



(C3) 
(C4) 



(G K A {\ 



G 



K 



Kj 
Bl 
Kj 
A2 



\G K B l 





+ |sj|e|7i 




e\ 


vW 2 - 


4 s jI £ I7i 




e 




+ 5«i|e|7i 


7(«j- 


-ky) 







\ , 

(C5) 

These evanescent modes exist in the region \yky\ > 
y/e 2 — Sj-|e|7i. Therefore, there are no traveling modes 
in both monolayer and bilayer graphenes in the region 

| 7 fc„| > V^+Wi- (C6) 

Note that G K 2 is the same as Eq. (fl"4)) . 

Edge states localized near the boundary (x 
the wave function 



G(r) = 



E 

v=K,K' 



a v G v _e 



- Kx+ik y y 



0) have 

(x < 0); 



E 

v=K,K' 3=1,2 



VJ G v +3 e-^ x+lk y y (x > 0), 



with appropriate coefficients a„ and /3„j 
for boundary ZZ1, we have 



(C7) 

More explicitly, 



/~lVl 
U A1 


U A1 


G V A \ 


/ Pv\ \ 








■nv2 


iG B 


Pv2 


= 0, 


(C8) 


rivl 
U B2 


r<v2 

U B2 


/ 


V -a v ) 







where we have multiplied imaginary unit i in such a way 
that the coefficient matrix becomes real. For boundary 
ZZ2, we have 



U A1 


U A1 


G\\ 


/ Pv\ \ 






;r<vl 
lU Bl 


■nv2 

lKj Bl 


iG% 


Pv2 


= 0, 


(C9) 


lU A2 


1< G V A2 


/ 


V —a v ) 







For AC1, we have 
1K1 



I U A1 

;nKl 
iu B1 

U B2 




A2 



U A1 



iG 



K2 
Bl 



r<K2 

U B2 






iG% 







G 
G 
iG 
iG 






K'l 
B2 
K'l 
Al 
K'l 
Bl 
K'l 
A2 






K'2 
B2 
K'2 
Al 
K'2 
Bl 

i G A 2 



G 
G 

iG 



G 
iG 



A" 
A 
A" 



J 



(3K2 

-a K 

P K'l 
Pk>2 

\-OL K > J 



12 



= 0. (CIO) 

The determinant of the coefficient matrix remains 
nonzero in the energy range satisfying Eq. (|C6[) and van- 
ishes at e = 0. Therefore, edge states can be present only 
at £ = 0. 

Let us consider the special case e = +0 or e = — 0. 
In the monolayer region, we have k = —k y for k y > 0, 
K = k y for k y < 0, and therefore the evanescent mode 
becomes 



G* = 



(k y > 0); 

(ky < 0), 



(Cll) 



where we have multiplied an appropriate phase factor. 
The wave function for the K' point is obtained by replac- 
ing k y with — k y . 

In the bilayer region, on the other hand, we have 
Kj = \k y \ for both j = 1 and 2, and consequently 
becomes the same between j = 1 and 2. In order to 
obtain two independent evanescent modes we expand 
G&(:c)=G£ 



- r ' K e KjX in terms of 5 = \e\-yi/(jk y )^ 



Gf + (x) = Gf + ( 



5Gf + (x 



0(6 2 



with 



G*Ax) 



S 3 Se 2 




1 + Sy 



-\ky\x 



(C13) 



and 



g£(x)( 



-SjS e 




l+£j, 
2 

2 

,7^1/1 



\ 



(l + 2\k v \x)e-^ x 



\SjS e 



7i 
7i 



V 



-|fcy|s 



(C14) 



/ 



Then, two independent modes can be written as G{^ (x) 
and Gf (x) with Gf (x) = GZ^xjW and 



Gf(x) 



£ Sj -G«»m (fc,>0); 

3=1,2 

£ G^-^W (fc,<0). 

3=1,2 



(C15) 



Therefore, we have for k y > 



GfOr) 




(C16) 



(C17) 



\ ^-\ky\xj 



and for fc y < 



Gf(z) 



G 2 K (x) 




(C18) 



(C19) 



(C12) The wave functions for the K' point are again obtained by 
replacing k y with —k y , i.e., G^ (x; k y ) = G^ (x; —ky). 
Therefore, we have for ky > 



ZZ1 



ZZ2 : 



AC1 : 



K : olk — Pki = Pk2 = 0, 
K' : a K > = 0, Pk'i + \Pk>2 

K : olk — Pk2 = 0, 
K' : a K > = Pk'I = Pk>2 = 0, 

OLK = pKl = PK2 = 0, 
OLK' = PK'l = Pk'2 = 0, 



0. 



and for k y < 



ZZ1 



ZZ2 



AC1 



K : a K =0, (3ki + ±(3k2 = 
K' : a K > = Pk'1 = Pk>2 = 0, 

K : a K = Pki = Pk2 = 0, 
K' : a K > = Pk>2 = 0, 

OK = Pk\ = PK2 = 0, 
OLK' = Pl<'l = Pk'2 = 0. 



(C20) 
(C21) 
(C22) 

(C23) 
(C24) 
(C25) 



There is no edge state in the armchair boundary. 

In the case of boundary ZZ1, we have a single edge 
state at the K' point for k y > and one at the K point for 
ky < 0. The wave function of these states is completely 
localized in the bilayer region and is given by 



G K (x) = 



( \K\ X \ 
o 

7i 

V o / 



(k y < 0). 



(C26) 
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The wave function for the K' point (k y > 0) is also given 
by the above equation. 

In the case of boundary ZZ2, on the other hand, we 
have a single edge state at the K point for k y > and 
one at the K' point for k y < 0. The wave function is 
completely localized in the bilayer region and is given by 



G K (x) 



{ky > 0). 



(C27) 



The wave function for the K' point (k y < 0) is again given 
by the same expression. These results are summarized in 
Table H 

For j — 2 and at |e| = 71, we have K2 — \k y \, giving 
Gf| = for k y < and Gf ' 2 2 = for k y > in Eq. 
(|C5[) . Other elements of and G j all remain nonzero. 
For the boundary ZZ1, therefore, the boundary condi- 



tion Fb2(0, y) = is satisfied and traveling modes in the 
monolayer can be connected only to the evanescent mode. 
It is easy to show that this evanescent mode cannot be 
connected to the evanescent mode in the monolayer and 
therefore cannot form a pure edge state. However, we 
have perfect reflection at \e\ = 71 for —71 < jk y < at the 
K point and for 0<jk y <ji at the K' point. 

This perfect reflection is closely related to the vanish- 
ing transmission probability at e = 71 for ZZ1 shown in 
Fig. [H In fact, when the Fermi level lies at 71 under the 
condition that the electron density is the same between 
the monolayer and bilayer graphenes, i.e., k — v / 27i/7j 
the reflection coefficient for wave incident from the mono- 
layer side is calculated as 



cos6> - i(l-2V2 sine ,„ . 

r KK = =r (6 > 0). 

cos 6 + t(l-2V2) sinS 



(C28) 
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